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Abstract. Interpretation of D-meson decay-dynamics has revealed itself to be strongly dependent 
' on our understanding of the light-meson sector. The interplay becomes particularly evident in Dalitz 

plot analyses to study physics within and beyond the Standard Model. Experience and results from 
FOCUS are presented and discussed. A brief update of the pentaquark search in the experiment is 
also reported. 
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INTRODUCTION 

Dalitz analyses are largely applied in modern high-energy experiments to study Heavy 
Flavor hadronic decays, but also to perform precise measurements of the CKM matrix 
elements and search for new physics. Paradigmatic examples are B — > pn and B — > 
! £)(*)£"(*) for the extraction of the a and 7 angles of the Unitarity Triangle. The extraction 

of a in B — ► pn requires filtering the desired intermediate states among all the possible 
combinations, e.g. on, /o(980)7T etc. The extraction of 7 in B — ► D^K^ 
requires, in turn, modeling the D amplitudes. The %% and Kit S-wave are characterized 
by broad, overlapping states: unitarity is not explicitly guaranteed by a simple sum 
of Breit-Wigner functions. In addition, independently of the nature of the o, it is 
not a simple Breit-Wigner. The /o(980) is a Flatte-like function, and its lineshape 
parametrization needs precise determination of KK and %% couplings. Recent analyses 
of CP violation in the B — ► DK channel from the beauty factories needed two ad hoc 
resonances to reproduce the excess of events in the tzk spectrum, one at the low-mass 
threshold, the other at 1 . 1 GeV 2 [Ql 0] • This procedure of "effectively" fitting data invites 
a word of caution on estimating the systematics of these measurements. A question then 
naturally arises: in the era of precise measurements, do we know sufficiently well how 
to deal with strong-dynamics effects in the analyses? We have faced parametrization 
problems in FOCUS and learnt that many difficulties are already known and studied in 
different fields, such as nuclear and intermediate-energy physics, where broad, multi- 
channel, overlapping resonances are treated in the K-matrix formalism [Hi 3, B] . The 
effort we have made consisted mainly in building a bridge of knowledge and language 
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to reach the high-energy community; our pioneering work in the charm sector might 
inspire future accurate studies in the beauty sector. FOCUS Dalitz plot analyses of the 
D+,D S — > Ji + n~ii + and of the D + — > K~iz + ii + will be discussed. 

The collaboration has also taken a complementary non-parametric approach to mea- 
suring the K tz + amplitude in the D + — > K K + n + decay using a projective weighting 
technique. Results will be presented. 



THE D + AND ZX 



7T+ AMPLITUDE ANALYSIS 



The FOCUS collaboration has implemented the K-matrix approach in the D s and D + — > 
Tt + K~K + analyses. Results and details can be found inJ6Q. It was the first application 
of this formalism in the charm sector. In this model pO, the production process, i.e, 
the D decay, can be viewed as consisting of an initial preparation of states, described 
by the F '-vector, which then propagates according to (I — iKp)~ l into the final one. 
The K-matrix here is the scattering matrix and is used as fixed input in our analysis. 
Its form was inferred by the global fit to a rich set of data performed in QVD - It is 
interesting to note that this formalism, beside restoring the proper dynamical features 
of the resonances, allows for the inclusion in D decays of the knowledge coming from 
scattering experiments, i.e, an enormous amount of results and science. No re-tuning of 
the K-matrix parameters was needed. The confidence levels of the final fits are 3.0% 
and 7.7 % for the D s and D + respectively. The results were extremely encouraging 
since the same K-matrix description gave a coherent picture of both two-body scattering 
measurements in light-quark experiments as well as charm-meson decay. This result was 
not obvious beforehand. Furthermore, the same model was able to reproduce features of 
the D + — ► n + n~n + Dalitz plot, shown in figfU that would otherwise require an ad hoc 
o resonance. The better treatment of the 5-wave contribution provided by the K-matrix 
model was able reproduce the low-mass k + k~ structure of the D + Dalitz plot. This 
suggests that any a-like object in the D decay should be consistent with the same a-like 
object measured in n + n~ scattering. 




Low m 2 / n - mass projection (GeV 2 /c 4 ) High 



mass projection (GeV 2 /c 4 ) Low m 2 * - mass projection (GeV 2 /c 4 ) High m 2 K + K - mass projection (GeV 2 /c 4 ) 



FIGURE 1. FOCUS Dalitz -plot projections for D s and D + to three pions with fit results superimposed. 
The background shape under the signal is also shown. 

Further considerations and conclusions from the FOCUS three-pion analysis were 
limited by the sample statistics, i.e. 1475 ±50 and 1527 ±51 events for D s and D + 
respectively. 

We considered imperative to test the formalism at higher statistics. This was accom- 
plished by the D + — > K~n + n + analysis. 



THE D+ -> K~7i + 7i + AMPLITUDE ANALYSIS 



The recent FOCUS study of the D + — > K~k + k + channel uses 53653 Dalitz-plot events 
with a signal fraction of ~ 97%, and represents the highest statistics, most complete 
Dalitz plot published analysis for this channel. Invariant mass and Dalitz plots are shown 
in figEl Details of the analysis may be found in dsj]. 




FIGURE 2. The D + — > K~n + n + Dalitz plot (left) and mass distribution (right). Signal and sideband 
regions are indicated: sidebands are at ±(6-8) a from the peak. 

An additional complication in the Kk system comes from the presence in the S- 
wave of the two isospin states, 1=1/2 and 1 = 3/2. Although only the / = 1/2 is 
dominated by resonances, both isospin components are involved in the decay of the D + 
meson into K~k + k + . A model for the decay amplitudes of the two isospin states can be 
constructed from the 2x2 K-matrix describing the 1=1/2 S- wave scattering in {Kn)\ 
and (Kr\') 2 (with the subscripts 1 and 2, respectively, labelling these two channels), and 
the single-channel K-matrix describing the 1 = 3/2 K K + — > K K + scattering. The K- 
matrix form we use as input describes the 5-wave K K + — > K K + scattering from the 
LASS experiment ||9|l for energy above 825 MeV and K K — > K K scattering from 
Estabrooks et al. Q 1QQ. The K-matrix form follows the extrapolation down to the Kk 
threshold for both 1=1/2 and 1 = 3/2 S-wave components by the dispersive analysis 
by Biittiker et al. [[Till , consistent with Chiral Perturbation Theory |fl2ll . The total D- 
decay amplitude can be written as 

^ = (F 1/2 ) 1 (s) +F 3/2 (s) +Y< a J e ' Sj B(abc\r), (1) 

j 

where s = M 2 (Kk) , (F { j 2 ) x and F 3 / 2 represent the / = 1 /2 and 1 = 3/2 decay amplitudes 
in the Kk channel, j runs over vector and spin-2 tensor resonances and B(abc\r) are 
Breit-Wigner forms. The J > resonances should, in principle, be treated in the same 
K-matrix formalism. However, the contribution from the vector wave comes mainly 
from the ^*(892) state, which is well separated from the higher mass ^f*(1410) and 
K* (1680) , and the contribution from the spin-2 wave comes from K^ (1430) alone. Their 



contributions are limited to small percentages, and, as a first approximation, they can be 
reasonably described by a simple sum of Breit-Wigners. More precise results would 
require a better treatment of the overlapping .ST* (1410) and ^* (1680) resonances as 
well. In accordance with SU(3) expectations, the coupling of the K% system to Kr\ is 
supposed to be suppressed. Indeed we find little evidence that it is required. Thus the 
F\j2 form for the K% channel is 

(^1/2)1 = {I-iK 1/2 p)$(P l/2 )j, (2) 

where / is the identity matrix, K± / 2 is the K-matrix for the 1=1/2 5-wave scattering in 
Kn and Kr]', p is the corresponding phase-space matrix for the two channels 140 and 
(P\/ 2 )j is the production vector in the channel j. 
The form for F 3 j 2 is 

F 3/2 = (I-iK 3/2 p)- l P 3/2 , (3) 

where K 3 / 2 is the single-channel scalar function describing the 1 = 3/2 K it + — > K Tt + 
scattering, and P3/2 is the production function into K%. 

The P-vectors are in general complex reflecting the fact that the initial coupling 

D+ -> [K- 

~ft + )ftf pec tator need not be real. Their functional forms are: 

(Pi ,2) 1 = + (cio + c u ? + c n ?)e'* (4) 

(^1/2)2 = + (c 2 o + C2is + c 22 f)e^ (5) 

P 3 /2 = (c30 + c 3 is + c 32 ^)e^. (6) 

/3e' e is the complex coupling to the pole in the 'initial' production process, gi and 
g 2 , s\ and s 2 are the K-matrix couplings and poles. The polynomials are expanded 
about T= s — s c , with s c = 2 GeV 2 corresponding to the center of the Dalitz plot. The 
polynomial terms in each channel are chosen to have a common phase ji to limit the 
number of free parameters in the fit and avoid uncontrolled interference among the 
physical background terms. Coefficients and phases of the P-vectors are the only free 
parameters of the fit determining the scalar components. Kn scattering determines the 
parameters of the K-matrix elements and these are fixed inputs to this D-decay analysis. 
Free parameters for vectors and tensors are amplitudes and phases (a,; and Table Q] 
reports our K-matrix fit results. It shows that quadratic terms in (Pi/ 2 )i are significant 
in fitting data, while in both (^1/2)2 an d P 3 / 2 constants are sufficient. The J > states 
required by the fit are listed in table |2l 

The S- wave component accounts for the dominant portion of the decay (83.23 ± 
1.50)%. A significant fraction, 13.61 ±0.98%, comes, as expected, from #*(892); 
smaller contributions come from two vectors ^*(1410) and ^*(1680) and from the 
tensor £|(1430). It is conventional to quote fit fractions for each component and this is 
what we do. Care should be taken in interpreting some of these since strong interference 
can occur. This is particularly apparent between contributions in the same-spin partial 



TABLE 1. 5-wave parameters from the K-matrix fit to the FOCUS 
D + — > K~n + n + data. The first error is statistic, the second error is 
systematic from the experiment, and the third is systematic induced by 
model input parameters for higher resonances. Coefficients are for the 
unnormalized S- wave. 

coefficient phase (deg) 

j3 = 3.389 ±0.152 ±0.002 ±0.068 9 =286±4±0.3±3.0 

ci = 1.655 ±0.156±0.010±0.101 J\ = 304±6±0.4±5.8 
a i = 0.780 ± 0.096 ± 0.003 ± 0.090 
ci2 = -0.954 ± 0.058 ± 0.0015 ± 0.025 

c 20 = 17.182±1.036±0.023±0.362 ^ = 126 ±3 ±0.1 ± 1.2 

c 30 = 0.7 34 ±0.080 ± 0.005 ±0.030 ^ = 211 ± 10 ±0.7 ±7.8 



Total S-wave fit fraction = 83.23 ± 1 .50 ± 0.04 ± 0.07 % 
Isospin 1/2 fraction = 207.25 ±25.45 ± 1.81 ± 12.23 % 
Isospin 3/2 fraction = 40.50 ±9.63 ±0.55 ±3.15 % 



TABLE 2. Fit fractions, phases, and coefficients for the J > compo- 
nents from the K-matrix fit to the FOCUS D + — > K~n + n + data. The first 
error is statistic, the second error is systematic from the experiment, and 
the third error is systematic induced by model input parameters for higher 
resonances. 



component fit fraction (%) phase 8j (deg) coefficient 



K*(892)k + 


13.61 ±0.98 


(fixed) 


1 (fixed) 




±0.01 ±0.30 






K*(m0)n + 


1.90 ±0.63 


1±7 


0.373 ±0.067 




±0.009 ±0.43 


±0.1±6 


±0.009 ±0.047 


/q(1430)7z; + 


0.39 ±0.09 


296 ±7 


0.169±0.017 


±0.004 ±0.05 


± 0.3 ± 1 


±0.010±0.012 


K*{Ul0)n + 


0.48 ±0.21 


293 ±17 


0.188±0.041 




±0.012±0.17 


±0.4±7 


±0.002 ±0.030 



wave. While the total 5-wave fraction is a sensitive measure of its contribution to the 
Dalitz plot, the separate fit fractions for / = 1/2 and 1 = 3/2 must be treated with care. 
The broad 1=1/2 S-wave component inevitably interferes strongly with the slowly 
varying 1 = 3/2 S- wave, as seen for instance in jl3\\ . Fit results on the projections are 
shown in fig. [31 

The fit % 2 I& .o.f is 1.27 corresponding to a confidence level of 1.2%. Our adaptive 
binning scheme is shown is fig [4] If the 1 = 3/2 component is removed from the fit, the 
^ 2 /d.o.f worsens to 1.54, corresponding to a confidence level of 10~ 5 . 

These results can be compared with those obtained in the effective isobar model, 
consisting in a sum of Breit Wigners, which can serve as the standard for fit quality. 
Two ad hoc scalar resonances are required, of mass 856 ± 17 and 1461 ±4 and width 
464 ± 28 and 177 ± 8 MeV/c 2 respectively to reproduce the data and reach a ^ 2 /d.o.f 
is 1.17, corresponding to a C.L of 6.8%. A detailed discussion of the results and the 
systematics can be found in A feature of the K-matrix amplitude analysis is that it 
allows an indirect phase measurement of the separate isospin components: it is this phase 




m 2 K - K + mass projection (GeV 2 /c 4 ) m 2 K - B + mass projection (GeV 2 /c 4 ) 

FIGURE 3. The Dalitz plot projections with the K-matrix fit superimposed. The background shape 
under the signal is also shown. 




FIGURE 4. The adaptive binning scheme. 



variation with isospin 7=1/2 that should be compared with the same 7=1/2 LASS 
phase, extrapolated from 825 GeV down to threshold according to Chiral Perturbation 
Theory. This is done in the right plot of fig. [51 In this model |5j] the P-vector allows for 
a phase variation accounting for the interaction with the third particle in the process of 
resonance formation. It so happens that the Dalitz fit gives a nearly constant production 
phase. The two phases in fig.[5h>) have the same behaviour up to ~ 1.1 GeV. However, 
approaching Kr\ ' threshold, effects of inelasticity and differing final state interactions 
start to appear. The difference between the phases in fig. [5k) is due to the 7 = 3/2 
component. 




FIGURE 5. Comparison between the LASS 1=1/2 phase + ChPT (continous line) and the F-vector 
phases (with ±1 a statical error bars); a) total F-vector phase; b) / = 1/2 F-vector phase. The vertical 
dashed line shows the location of the Kr]' . 



These results are consistent with Kn scattering data, and consequently with Watson's 
theorem predictions for two-body K% interactions in the low K% mass region, up to ~ 1.1 
GeV, where elastic processes dominate. This means that possible three-body interaction 
effects, not accounted for in the K-matrix parametrization, play a marginal role. The K- 
matrix form used in this analysis generates the S-matrix pole E = M — iF/2 = 1.408 — 
/0.110 GeV. Any more distant pole than ^(1430) is not reliably determined as this 
simple K-matrix expression does not have the required analyticity properties. However, 
our K-matrix representation fits along the real energy axis inputs on scattering data and 
Chiral Perturbation Theory in close agreement with those used in ||14I1 . which locates 
the K with a mass of (658 ± 13) MeV and a width of (557 ± 24) MeV by careful 
continuation. These pole parameters are quite different from those implied by the simple 
isobar fits. We have thus shown that whatever K is revealed by our D + — > K~n + n + 
results, it is the same as that found in scattering data. 



A NON-PARAMETRIC APPROACH TO DETERMINE THE K~ti + 
AMPLITUDE IN D+ -> K~K + 7i + DECAY 



While making the effort of refining the amplitude formalism, FOCUS identified the 
D + — > K K + n + as an ideal case to apply the projective weighting technique developed 
in the semi-leptonic sector ft 1 311 to the hadronic decays, with no need to assume specific 
Breit-Wigner resonances, forms for mass dependent widths, hadronic form factors or 
Zemach momentum factors. Details can be found in ||16tl . The old E687 Dalitz plot anal- 
ysis IU7I1 concluded that the observed D + — ► K~K + iz + Dalitz plot could be adequately 
described by just three resonant contributions: (f) n + ,K + K* (892) and K + K*o(1430). Al- 



though (j)7C + i s an important contribution, the is a very narrow resonance that can 
be substantially removed through a cut on m K + K -, i.e m K + K - > 1050 MeV/c 2 . Since 
there is no overlap of the band with the K* and most of the kinematically allowed 
^*o(1430) region, there is a relatively small loss of information from the anti-0 cut; of 
course careful systematic evaluation for residual K + K contributions and bias are per- 
formed. In the absence of the K K + resonances, we can write the decay amplitude in 
terms of m K - K + = m and the helicity decay angle 6 Thus 

s,p.d... 

A= £ A,(m)4)(cos0), (7) 
/ 

where d l 00 (cos6) are the Wigner d-matrices describing the amplitude for a K 7t + system 
of angular momentum / to simultaneously have angular momentum along its helicity 
axis and the K n + decay axis. This technique is an intrinsically one-dimensional anal- 
ysis. The decay intensity assuming, for simplicity, that only S and P-waves are present, 
is 



\A\ 2 = \S(m) +P(m) cos 9\ 2 = \S(m) | 2 + 2Re{S*(m)P(m)}cos 9 + \P(m) | 2 cos 2 0, (8) 
The approach is to divide cos into twenty evenly spaced angular bins. Let 

~fi = ( i n h i n 2 - i n 2 o) (9) 

be a vector whose 20 components give the population in data for each of the 20 cos 
bins. Here ' specifies the i th m K - K + bin. Our goal is to represent the D vector in eq. |9] 
as a sum over the expected populations for each of the three partial waves. For this 
simplified case there are three such vectors computed for each m K - K + bin, {'Jft a } = 
( l nt ss ,'Jttsp, l ntpp)- 

Each l ~nt a is generated using a phase- space and full detector simulation for D + — > 
K K + lt + decay with one amplitude turned on, and all the others turned off. We then 
use a weighting technique to fit the bin populations in the data to the form 'D = 
Fss(mi) l ~trf ss + Fsp(mj)'nt sp + Fpp(mi) l lft pp. When including the D-wave as well, the 
results appear just as five weighted histograms in the m K - K + mass, as in fig. [6l for the 
five independent amplitude contributions. 

The curves in fig.[6]are the model used in E687 but with a ^*o( 1430) ad hoc arranged 
to fit the data, i.e represented as a Breit Wigner with a pole at mo=1412 MeV/c 2 
and a width of T = 500 MeV/c 2 , not consistent with the standard PDG i^*o(1430) 
parametrization used by E687. This analysis reveals once more, how subtle the inclusion 
of the broad S-wave resonances in charm Dalitz analysis can be. Although the D + — > 
K K + tz + is an ideal case, it might be possible to extend the analysis to the D s — ► 
K K + n + decay, as well as D° — ► K + K K° and hadronic four body decays such as 
D° -> K K+Ti+Ti -> tyn+ir. 
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FIGURE 6. This figure compares the five projected amplitudes obtained according to their angular 
dependence (error bars) with the E687 model properly corrected, as explained in the text, to match the 
data (red curves). The plots are: a) S 2 direct term, b) 2 S xP interference term, c) P 2 direct term, d)2 PxD 
interference term and e) D 2 direct term. 



SEARCH FOR PENTAQUARK CANDIDATES 



The FOCUS collaboration searched for the charmed <d Q c (cuudd) pentaquark candidate in 



the decay modes 0^ — > D* p and ®®. — > D~p II 1 311 . No evidence for a pentaquark at 3.1 
GeV/c 2 or at any mass less than 4 GeV/c 2 was observed. More recently the search was 
extended to two other candidates: & + (suudd) — > pK® Iil9ll and § (lS60)(ssddu) — > 
Z~n~ |[2~oll . Having found no evidence, limits were calculated. The + production cross 
section was normalized to E*(1385) ± and K*(S92) + because the reconstructed decay 
modes of the particles E*(1385) ± — >■ A°^ ± and K*(&92) + -> K^7t + are very similar, 
in terms of topology and energy release, to the signal. The 95 % C.L upper limits of 

° {@ l\™^f s) < 0.00012 (0.00029) and ^"Sss^ < °- 0042 (°- 0099 ) were 
estimated for a natural width of and 15 MeV/c 2 in the good acceptance region of 
the detector , i.e. for parent particles with momenta above 25 GeV/c. Analogously 
the upper limit was calculated for the Sj~(0 (I860)) candidate with respect to the 

S*(1530)° -> Z-n + obtaining g(Ss SSp-^ < °- 007 (°- 019 ) for a natural width 
of 0(1 5) MeV/c 2 . 



CONCLUSIONS 



Dalitz-plot analysis represents a unique, powerful and promising tool for physics studies 
within and beyond the Standard Model; however to perform such sophisticate analyses, 
we need to model the strong interaction effects. FOCUS has performed pilot studies in 
the charm sector through the K-matrix formalism and has started an effort to identify 
channels where non-parametric approaches can be undertaken. What has been learnt 
from charm will be beneficial for future accurate beauty measurements. 
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